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ABSTRACT  •+ 


An  experimental  and  analytical  program  was  conducted  to  develop  a 
methodology  to  predict  total  recession,  including  the  contribution  by 
particulate  impact  erosion,  in  solid  propellant  rocket  nozzles.  The  program 
was  divided  into  two  phases.— i  Under  Phase  d,'  a  methodology  was  developed  for 
predicting  recession  due  to  particle  impact  erosion  in  the  expansion  (exit 
cone)  region  of  nozzles. s  This  phase  was  performed  by  the  Aerojet  Strategic 
Propulsion  Company  (ASPC)  under  a  subcontract  to  the  Aerotherm  Division  of 
Acurex  Corporation.  Under  Phase  Iir  a  methodology  was  developed  for  predicting 
recession  due  to  erosion  coupled  with  ablation  in  the  entrance  and  throat 
regions  of  nozzles^ — This  phase  was  performed  by  Aerotherm. 

For  the  supersonic  expansion  region  methodology  (Phase  I),  the  TD2P 
computer  code,  contained  within  the  Interim  Solid  Performance  Program  (ISPP), 
was  chosen  for  predicting  the  two-phase  flow.  The  TD2P  code  is  used  to 
predict  the  location  of  particle  impingement  on  the  nozzle  exit  cone  wall 
and  the  impact  conditions.  Also,  it  was  used  in  conjunction  with  erosion  test 
data  to  generate  mechanical  erosion  models.  Data  for  erosion  models  were 
obtained  from  one  series  of  tests  in  the  Dust  Erosion  Tunnel  (DET)  at  the  Arnold 
Engineering  Development  Center  (AEDC)  and  a  second  series  of  tests  in  small 
rocket  motors  at  ASPC.  Two  independent  mechanical  erosion  models  were  developed. 
The  first  was  derived  using  the  "G-law"  approach,  while  the  second  was  formulated- 
to  represent  the  physical  phenomena  occurring  during  the  erosion  process. 

For  the  subsonic  entrance  and  throat  region  methodology  (Phase  II) ,  a 
variety  of  tests  were  used  in  modeling  particle  impact  erosion  —  particle 
(droplet)  breakup ;  aluminum  (Al)  and  aluminum  oxide  (Al^O^)  chemical  reactions 
with  carbon;  and  subsonic  erosion  tests  using  Al  and  Al^O^  impacting  a  carbon 
surface.  The  subsonic  erosion  was  modeled  using  the  "G-law"  approach.  Two 
computer  codes  were  developed:  (1)  the  Chamber  Flowfield  Code  (CFC)  for 
predicting  particle  impact  conditions  at  the  nozzle  wall,  and  (2)  the  Charring 
Material  Ablation-Erosion  (CMAE)  code  for  computing  the  contributions  of  erosion 
and  ablation  to  recession. 

Both  the  Phase  I  and  Phase  II  recession  prediction  methodologies  were 
demonstrated  by  comparing  predictions  to  data  from  rocket  motor  test  firings 
having  known  or  suspected  particle  impact  erosion.  Program  results  are  reported 
in  four  volumes,  as  follows: 

•  Volume  I:  Recession  Prediction  Methodology  for  Rocket  Nozzle 

Entrance  and  Throat  Regions 

•  Volume  II:  Erosion  Model  Development  (Expansion  Region) 

•  Volume  III:  User's  Manual  —  Chamber  Flowfield  Code  (CFC)  and 

Charring  Material  Ablation-Erosion  (CMAE)  Code 

•  Volume  IV:  User's  Manual  —  Erosion  Prediction  Procedure  for 

Rocket  Nozzle  Expansion  Region 


SUMMARY 
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An  experimental  and  analytical  program  mat  conducted  to  develop  a 
methodology  to  predict  total  recession,  including  the  contribution  by 
particulate  impact  erosion,  in  solid  propellant  rochet  nozzles.  The  program 
was  divided  into  two  phases.  Under  Phase  I,  a  methodology  was  developed  for 
predicting  recession  due  to  particle  impact  erosion  in  the  expansion  (exit 
cone)  region  of  nozzles.  This  effort  was  performed  by  the  Aerojet  Strategic 
Propulsion  Company  (ASPC)  under  a  subcontract  to  the  Aerotherm  Division  of 
Acurex  Corporation.  Under  Phase  II,  a  methodology  was  developed  for  predicting 
recession  due  to  erosion  coupled  with  ablation  in  the  entrance  and  throat 
regions  of  nozzles.  This  effort  was  performed  by  Aerotherm.  Results  from  this 
program  are  reported  in  four  separate  volumes,  as  follows: 


Volume  I: 


.  Volume  II : 


.  Volume  III: 


Volume  IV: 


Recession  Prediction  Methodology  for  Rochet  Nozzle 
Entrance  and  Throat  Regions 

Erosion  Model  Development  (Expansion  Region) 

User's  Manual  —  Chamber  Flowfield  Code  (CFC)  and 
Charring  Material  Ablation-Erosion  (CMAE)  Code 

User's  Manual  —  Erosion  Prediction  Procedure  for 
Rochet  Nozzle  Expansion  Region 
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lht  remainder  of  this  summary  provides  a  brief  overview  of  the 
recession  prediction  methodologies  developed  for  the  nozzle  expansion  region 
(reported  in  Volume  II)  and  for  the  nozzle  entrance  and  throat  regions 
(reported  in  Volume  I). 


t 
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Nozzle  Expansion  Region 

A  procedure  was  developed  for  predicting  the  erosion  that  occurs  in  a 
solid  rochet  nozzle  exit  cone  when  aluminum  oxide  (AlgO^)  particles  in  the  *"low 
inpact  the  nozzle  wall  in  the  supersonic  region  downstreaa  of  the  nozzle 
throat.  This  procedure  consists  of  a  two-phase  flow  solution  and  two 
independent  nechanical  erosion  models.  These  models  calculate  the  nozzle  wall 
erosion  patterns  in  carbon-carbon  and  bulk  graphite  materials  as  a  function  of 
the  impinging  particle  flow  properties,  the  impact  angle,  and  tho  nozzle  wall 
material  properties. 
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An  empirical  'G-law'  erosion  model  was  developed  for  nozzle  expansion 
regions  from  Dust  Erosion  Tunnel  (DET)  test  data  and  small  solid  motor  firing 
data.  Specimens  of  carbon-carbon,  bulk  graphite,  and  tantalum/ tungsten  were 
exposed  to  particle-laden  test  gases.  A  debris  layer  shielding  erosion  model 
was  also  developed  for  the  expansion  region.  In  contrast  to  the  empirical 
G-law,  this  model  was  formulated  to  represent  the  fundamental  physics  of  the 
erosion  process. 

The  TD2P  two-phase  flow  code  contained  in  te  Interim  Solid  Performance 
Program  (ISPP)  was  modified  to  incorporate  the  above  two  erosion  models.  This 
version  of  TD2P  constitutes  the  methodology  for  predicting  particle  impact 
conditions  at  the  walls  of  nozzle  exit  cones  and  the  resulting  wall  erosion. 
Predictions-  were  made  using  both  erosion  models  for  various  solid  rocket  motor 
test  firings  where  erosion  occurred.  Predictions  using  the  G-law  model  were 
generally  acceptable  relative  to  test  results,  providing  the  particle  flow  and 
impact  conditions  were  within  the  range  covered  by  the  DET  and  small  motor 
experiments.  Predictions  using  the  debris  layer  model  wero  also  generally 
satisfactory,  but  not  as  accurate  as  G-law  predictions  when  particle  impact 
conditions  were  within  the  DET  and  small  motor  test  range.  However,  the  debris 
layer  model  was  the  most  accurate  for  conditions  outside  the  experimental 
range. 
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Entrance  and  Throat  Regions 


A  aethedology  was  developed  for  predicting  the  total  recession, 
including  particle  impact  erosion,  for  the  entrance  and  throat  regions  of  solid 
rocket  nozzles.  Major  program  efforts  included  the  development  of  a  motor 
chamber-nozzle  flowfield  code  and  a  model  for  describing  the  erosion  process  in 
the  subsonic  regions  of  the  nozzle* 

The  Fully  Coupled  Transonic  (FCT)  two-phase  flow  module  of  the  Solid 
Performance  Prediction  (SPP)  program  was  modified  to  compute  particle 
trajectory  crossovers  and  particle  impact  conditions  at  the  nozzle  wall.  Also, 
the  calculation  domain  was  extended  forward  to  account  for  the  motor  chamber 
and  propellant  grain,  including  mass  injection  from  the  grain.  However,  flow 
within  the  submerged  region  of  nozzles  remains  to  be  included.  This 
modification  is  referred  to  as  the  Chamber  Flowfield  Code  (CFC). 


Empirical  subsonic  erosion  models  were  developed  for  predicting 
recession  due  to  the  impact  of  aluminum  (Al)  and  Al2<>3  particles.  Models  were 
-  constructed  in  the  form  of  a  G-law  with  two  components:  mechanical  erosion  due 
to  mechanical  mass  removal  and  chemical  erosion  due  to  particle-surface 
chemical  reactions.  However,  the  models  are  based  on  bulk  graphite  data  and 
provide  only  rough  estimates  for  the  erosion  of  carbon-carbon  composites.  The 
Charring  Material  Ablation  (CXA)  code  was  modified  to  include  these  models  in 
predictions  of  recession  resulting  from  ablation  and  erosion.  This 
modification  is  designated  the  Charring  Material  Ablation-Erosion  (CMAE)  code. 


The  recession  prediction  methodology  for  subsonic  flow  regimes  was 
demonstrated  for  two  motor  firings  having  suspected  particle  erosion  in  the 
entrance  region  of  carbon-carbon  composite  nozzle*.  The  demonstration  results 
substantiated  the  prediction  of  the  onset  of  particle  impact  for  the  two 
motors.  However,  recession  predictions  using  the  bulk  graphite  erosion  model 
provided  only  a  rough  estimate  of  material  removed. 
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Although  not  currently  included  in  the  recessic^  prediction 
methodology,  particle  breakup  in  nozzle  entrance  and  throat  regions  was  modeled 
in  terms  of  Bond  number  and  nondimens ional  time.  This  model  was  based  upon 
cold-flow  data,  and  its  applicability  to  hot-firing  conditions  has  not  been 
demonstrated. 
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INTRODUCTION 


1.0 


This  Users  Manual  contains  the  information  needed  to  prepare  input 
for,  run,  and  interpret  results  from  two  separate  computer  codes  which 
calculate  erosion  due  to  particle  impact.  The  two  codes  are: 


1.  A  modified  version  of  the  Interim  Solid  Performance  Program 
(ISPP),  which  computes  the  two-phase  flow  in  a  rocket  nozzle 
and  thus  the  particle  impact  conditions,  if  any.  These 
conditions  are  then  used  to  calculate  the  resulting  nozzle 
wall  erosion. 

2.  A  stand-alone  code,  PIEROS,  in  which  the  particle  impact 
conditions  must  be  input  before  the  erosion  is  calculated. 

Both  codes  contain  two  distinct  particle  impact  erosion  models:  1)  an 
empirical  model  based  solely  on  experimental  data,  and  2)  a  debris-layer  model 
that  represents  the  actual  physical  processes.  An  option  is  available  for 
saving  the  particle  impact  conditions  calculated  during  an  ISPP  run  to  be  used 
as  subsequent  input  to  the  stand-alone  PIEROS  code.  The  various  options 
available  for  using  these  codes  are  represented  in  a  flow  diagram  presented  in 
Figure  C-l. 


When  using  the  ISPP  to  calculate  the  particle  flow  field  in  a  nozzle, 
the  erosion  resulting  from  any  particle  impact  at  the  nozzle  wall  is 
automatically  calculated  with  the  empirical  erosion  model  and  the  results 
printed  at  each  impact  location.  Since  the  debris-layer  erosion  model  requires 
substantially  more  computer  time  (usually  3  to  5  minutes),  its  use  during  the 
ISPP  run  is  left  as  an  option.  To  allow  subsequent  use  of  the  debris-layer 
model  (and/or  the  empirical  model  with  different  impact  parameters),  the 
particle  impact  conditions  calculated  with  the  ISPP  can  be  saved  on  a  unit 
specified  by  the  input  parameter  ISAVE. 


The  stand-alone  particle  impact  erosion  code,  PIEROS,  can  be  used 
independently  from  the  ISPP  code  by  having  the  particle  impact  conditions  input 
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A.  PARTICLE  FLOW  FIELD  AND  PARTICLE  B.  PARTICLE  IMPACT  EROSION 

IMPACT  EROSION  CALCULATION  CALCULATION  ONLY 


Figure  C-l.  Particle  Impace  Erosion  Calculation  Options 


directly,  or  it  can  use  the  conditions  saved  during  a  previous  ISPP  run.  As 
vith  the  ISFP,  particle  impact  erosion  in  PIEROS  is  calculated  automatically 
with  the  empirical  erosion  model,  with  use  of  the  debris-layer  model  optional. 

1.1  INTERFACE  WITH  ISPP 

ISPP  is  a  collection  of  computer  codes  designed  to  calculate  the 
delivered  specific  impulse  of  a  solid  rocket  motor.  The  overall  flow  diagram 
for  the  ISPP  is  shown  in  Figure  C-2.  A  complete  ISPP  run  would  include  the 
following  calculations: 

1.  A  one-dimensional,  shifting-equilibrium  nozzle  expansion  for 
a  specified  propellant  chemistry  with  the  ODE  code, 

2.  A  solid  propellant  grain  design  and  corresponding  internal 
ballistics  as  a  function  of  time  with  the  BAL  code, 

3.  A  one-dimensional,  finite-rate  nozzle  expansion  with  the  ODK 
code, 

4.  A  two-dimensional,  two-phase  nozzle  expansion  with  the  TD2P 
code, 

5.  A  turbulent  boundary  layer  solution  along  the  nozzle  wall  with 
the  TBL  code. 

However,  for  particle  impact  erosion  purposes  only  the  TD2P  code  is 
required,  possibly  in  conjunction  with  the  ODE  code  to  provide  the  flow 
properties  for  TD2P.  Fortunately,  ISPP  is  structured  to  allow  this  particular 
combination,  as  indicated  by  the  asterisks  in  Figure  C-2. 

It  is  not  intended  in  this  Users  Manual  to  include  the  details  of 
using  the  ISPP  codes,  which  are  covered  in  Reference  C-l.  However,  these  areas 
which  have  a  direct  bearing  on  the  particle  impact  erosion  calculations  will  be 
discussed.  It  is  assumed  tha?  the  user  already  has  the  ISPP  operational  on  a 
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DRIVER 


MAIN  PROCRAM 


INPUT  THERMOOY NAMIC  DATA  NOT  CONTAINED 
ON  TAPE  OR  PERMANENT  FILE 


INPUT  THERMODYNAMIC  DATA  FOR  TEMPERATURES 
BELOW  30Q°K  (IF  NEEDED) 


II  READ  SSPP  NAMELIST  (ODE.  BAL,  ODK,  TD2P.  TBL.  NOZZLE  GEOMETRY) 

2)  CALCULATE  NOZZLE  WALL  COORDINATES 

3)  READ  LINKAGE  DATA  FROM  PREVIOUS  RUNS  (IF  ANY) 

1)  REAOSOOE  NAMELIST  (PROPELLANT  CHEMICAL  COMPOSITION) 

2)  CALCULATE  ONE  DIMENSIONAL.  SHIFTING  EQUILIBRIUM  NOZZLE 
EXPANSION 

A)  THEORETICAL  SPECIFIC  IMPULSE 

B)  NOZZLE  FLOW  PROPERTIES 

II  REAOSBAL  NAMELIST  (GRAIN  DESIGN  AND  rfALLISTICS  INPUT) 

2)  CALCULATED  GRAIN  DESIGN  (IF  NEEDED) 

3)  INTERNAL  BALLISTICS  PARAMETERS  VS  TIME 

1)  READ  SODK  NAMELIST  (FiNITE  RATE  CHEMICAL  REACTION  DATA) 

2)  CALCULATE  ONE-OIMENSIONAL.  FINITE  RATE  NOZZLE  EXPANSION 


1)  READ  STD2  NAMELIST  (FLOW  PROPERTIES  (IF  NOT  AVAIL.  FROM 

ODE  OR  LINKAGE  DATA).  PARTICLE  DATA.  TRANSONIC  ANALYSIS  INPUT) 

21  CALCULATE  TWO-DIMENSIONAL,  TWOPHASE  NOZZLE  '’XPANSION 

II  READSTBL  NAMELIST  (FLOW PROPERTIES  (IF  NOT  AVAIL.  FROM 

ODE  OR  LINKAGE  DATA).  EDGE  CONDITIONS  (IF  NOT  AVAIL.  FROM  TD2P 
OR  LINKAGE  DATA) 

2)  CALCULATE  TURBULENT  BOUNDARY  LAYER  GROWTH  ALONG  NOZZLE 
WALL 


USED  IN  NORMAL  PARTICLE  IMPACT  EROSION  CALCULATION 
MAY  BE  USED  TO  CALCULATE  FLOW  PROPERTIES 


Figure  C-2.  ISPP  Flow  Diagram 


local  computer  system  and  that  the  modified  and  new  subroutines  needed  for  the 
particle  impact  erosion  calculations  can  simply  be  linked  to  the  existing  main 
program  and  associated  subroutines. 

As  shown  in  Figure  C-2,  a  typical  ISPP  calculation  for  particle 
impact  erosion  would  require  imput  from  the  $SPP  and  $TD2  namelists.  The 
nozzle  geometry  and  the  option  parameters  (ODE,  BAL,  ODK,  TD2P,  and  TEL)  are 
input  in  the  $SPP  namelist,  and  the  flow  properties,  particle  data,  and 
transonic  analysis  parameters  in  the  $TD2  namelist.  If  the  ODE  code  is  to  be 
used  for  calculating  the  flow  properties,  instead  of  inputting  them  directly,, 
then  the  necessary  data  is  input  in  the  $ODE  namelist,  preceded  by  a  list  of 
reactants  specifying  the  propellant  chemical  composition  (see  Reference  C-l  for 
details) . 

To  integrate  the  particle  impact  erosion  calculations  into  the  TD2P 
two-phase  flew  nozzle  solution,  four  existing  subroutines  were  modified  and 
four  new  subroutines  added.  A  flow  diagram  sh-dng  the  most  significant 
subroutines  used  in  this  solution  is  given  in  Figure  C-3. 

The  controlling  subroutine,  TD2P,  first  calls  INPTD2  to  read  input 
from  the  $TD2  namelist,  and  then  DIST  to  calculate  the  particle  size-varsus 
-mass  distribution.  Two  transonic  solutions  arc  available  for  the  two-phase 
flow  through  the  nozzle  inlet  and  throat,  the  Approximate  Transonic  Analysis 
(ATA)  or  the  Fully-Coupled  Transonic  (FCT).  The  choice  of  transonic  solution 
is  determined  by  the  parameters  ATA  and  FCT  from  the  $TD2  namelist.  The 
transonic  solution  terminates  just  downstream  of  the  throat,  generating  a  start 
line  across  the  nozzle  for  the  method-of-characteristics  (HOC)  solution  used 
for  the  supersonic  flow  in  the  nozzle  exit  cone. 

After  the  transonic  solution,  the  TD2P  subroutine  calls  N2MAIH,  which 
in  turn  calls  N3MAIN  to  integrate  the  thrust  across  the  startline,  and  CNTRL 
for  the  HOC  solution.  At  each  wall  solution  point  in  the  HOC  solution,  CNTRL 
calls  the  subroutine  THRST  to  integrate  the  thrust  along  the  wall  from  the 
previous  wall  solution  point,  and  the  subroutine  IMPNG  to  check  for  particle 
impingement  at  the  wall. 


To  include  the  particle  impact  erosion  calculations,  modifications 
were  made  to  the  TD2P,  DIST,  THRST,  and  IMPNG  subroutines,  and  the  ERODE 
subroutine,  including  the  entry  ERODEA,  was  added  along  with  supporting 
subroutines  DEBRIS,  LOOK,  and  KELOSS. 

To  allow  the  calculated  particle  impact  conditions  to  be  saved  on 
unit  ISAVE,  it  was  also  necessary  to  make  a  single  change  to  the  main  program 
DRIVER.  This  change  terminated  the  ISPP  run  at  the  conclusion  of  each  case, 
allowing  the  file  created  on  unit  ISAVE  to  be  closed. 

1.1.1  Modifications  to  Existing  Subroutines 

1.1. 1.1  DIST  Subroutine 

Originally  the  particle  size-vs-mass  distribution  could  be 
either  input  directly  or  calculated  by  the  DIST  subroutine,  with  the  mass  mean 
diameter  and  standard  deviation  either  input  or  calculated  internally.  The 
calculated  distribution,  however,  was  limited  to  equal  mass  increments;  e.g., 
if  5  particle  sizes  were  requested,  the  sizes  were  calculated  with  each 
representing  20  percent  of  the  particle  mass.  For  particle  impact  erosion,  it 
was  found  that  a  better  representation  of  the  smaller  particle  sizes  was 
needed,  which  could  only  be  obtained  by  requesting  a  large  number  of  particle 
sizes  with  a  corresponding  large  increase  in  computer  run  time.  Thus,  the  DIST 
subroutine  was  modified  so  that  unequal  mass  increments  could  be  input  and  the 
corresponding  particle  sizes  determined.  For  example,  6  particle  sizes 
representing  particle  mass  fractions  of  .01,  .02,  .04,  .08,  .15,  and  .70  can 
now  be  used. 

1.1. 1.2  IMPNG  and  THRST  Subroutines 

The  IMPNG  subroutine  was  modified  to  make  the  initial  call  to  the 
particle  impact  erosion  subroutine  ERODE  when  particle  impingement  is  first 
detected.  On  the  first  call  to  ERODE,  the  namelist  $ER0DIN  is  read  with  the 
data  required  for  the  erosion  calculations.  The  THRST  subroutine  was  modified 
to  call  ERODE  at  subsequent  wall  points  downstream  of  the  first  impact  point. 
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The  use  of  IMPNG  and  THKST  to  call  ERODE  is  repeated  for  each  particle  size 
that  impacts  the  nozzle  wall. 

1.1. 1.3  TD2P  Subroutine 

The  fourth  existing  subroutine  modified  wss  TD2P  to  print  a  summary 
of  the  calculated  particle  impact  erosion  at  the  conclusion  of  the  HOC 
solution,  and  to  call  the  entry  ERODEA  within  the  ERODE  subroutine  to  calculate 
the  erosion  with  the  debris-layer  model,  if  desired. 

1.1.2  Hew  Subroutines 


The  ERODE  subroutine,  new  to  1SPP,  reads  the  $ERODIN  namelist  and 
calculates  the  erosion  depth  at  each  point  along  the  wall  where  particle  impact 
occurs.  The  calculated  erosion  results  are  then  immediately  printed.  The 
particle  flow  conditions  at  each  impact  point  are  saved  for  subsequent  use  by 
the  debris-layer  model  and,  if  desired,  are  written  on  a  specified  output  file. 

If  the  debris-layer  model  is  to  be  used  for  calculating  erosion,  the 
entry  ERODEA  within  ERODE  is  called  by  TD2P  after  the  nozzle  flow  solution  has 
been  completed.  The  debris-layer  model  requires  an  axial  distribution  of 
particle  flow  properties  from  the  initial  impact  point  to  the  nozzle  exit. 
Three  subroutines,  also  new  to  ISPP,  are  called  from  ERODEA:  DEBRIS,  which 
calculates  flow  conditions  within  the  debris  layer,  LOOK,  an  interpolation 
routine,  and  KELOSS ,  which  calculates  the  kinetic  energy  lost  for  each 
collision  within  the  debris  layer. 

To  summarize,  any  current  ISPP  code  can  be  updated  to  include  the 
particle  impact  erosion  calculates  by  linking  the  updated  DRIVER  main  program 
and  the  following  modified  (m)  or  new  (n)  subroutines  to  the  existing  code: 
TD2P(m),  DIST(m) ,  THRST(m),  and  DCPNG(m),  ERODE(n) ,  DEBRIS(n),  LOOK(n),  and 
KELOSS(n). 
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2.0 


INPUT  DESCRIPTION 


2.1  ISPP  OPTION 


2.1.1  ISPP  Input 

2. 1.1.1  $SPP  Namelist 

No  changes  have  been  made  to  the  input  information  contained  vithin 
the  $SPP  namelist  (see  Reference  G-l).  To  run  only  the  TD2P  module: 

ODE  *  0. 

SAL  -  0. 

ODK  -  0. 

TBL  -  0. 

TD2P  -  1. 

If  the  ODE  is  to  be  used  for  calculating  the  flow  properties  for 

TD2P: 


ODE  «  1 

If  the  results  from  a  previous  ODE  calculation  are  to  be  used: 

ODE  ■  2 

The  only  other  input  data  required  from  the  $SPP  namelist  is  the 
nozzle  geometry  (pages  2-1  through  2-20  in  Reference  C-l).  If  plots  of  the 
method  of  characteristics  solution  along  with  the  gas  or  particle  streamlines 
are  desired,  the  plot  data  described  in  pages  2-25  through  2-30  in  Reference 
C-l  should  be  input.  The  particle  streamline  plots  can  be  useful  in  evaluating 
the  particle  impact  results. 
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2. 1.1. 2  $TD2  Namelist 


For  a  particle  impact  erosion  calculation,  the  $TD2  namelist  input 
described  on  pages  2-70  through  2-80  of  Reference  C-l  remains  essentially  the 
same.  The  only  change  is  the  additional  option  now  available  for  calculating 
the  particle  size-vs-mass  distribution.  This  option  is  invoked  by  inputting 
the  particle  mass  fractions  WPWT  (page  2-74)  without  the  corresponding  particle 
radii  R  (page  2-73).  Mote  that  the  number  of  WPWT  values  input  must  correspond 
to  MFG,  the  number  of  particle  size  groups,  and  that  the  summation  of  the  WPWT 
values  must  equal  1.0. 

2.1.2  Particle  Impact  Erosion  Input 

The  additional  input  required  for  particle  impact  erosion 
calculations  is  contained  in  the  namelist  $ER0DIN,  which  follows  the  $TD2 
namelist.  The  input  parameters  are  described  below: 

Parameter  Description  Onits  Default  Value 

$ER0DIN 

For  Debris-Laver  Erosion  Model: 


ST 

Yield  shear  stress  of  target  material 

psi 

2200. 

RHOTAR 

Target  material  density 

lbm/cu-ft 

120. 

FACTOR 

Coefficient  in  single  particle  erosion 
equation  (C  in  equation  4-b  in  report) 

none 

6.782 

EXPFAC 

Velocity  exponent  in  single  particle 
equation  (C2  in  equation  4-b  in  report) 

none 

2.0 

CD 

Particle  drag  coefficient  within  debris 
layer 

none 

5.0 

MART 

Maximum  number  of  finite  elements 

(calculation  terminates  when  exceeded) 

none 

1000 

IEMAX 

Maximum  number  of  consecutive  steps 

allowed  with  calculated  erosion  rate 

less  than  0.01  mil/sec 

none 

100 

1 


ERKIN 

Minimal  erosion  rate  (calcnlated  with 

eapirical  aodel)  that  anst  be  exceeded 

before  debris-layer  aodel  calcnlation 

begins 

ail/sec 

0. 

IPRINT 

Axial  distance  print  inereaent  in  debris- 

layer  aodel  (set  to  zero  to  delete  debris 

layer  print,  except  for  snaaary) 

inches 

0. 

For  Eapirical  Erosion  Model: 

T  -  OONSIE  (V_)  ■[,(  )*™»  (a  J6*’™  (d  yam [!-•»» (-EffW/P.V  > ] 

,  ^  P  P  P  P 

where  y  «  erosion  rate  in  ail/sec. 

V  «  psrticle  velocity  in  ft/sec, 

P  2 

ppVp  *  particle  aass  flux  in  lba/ft  ‘sec, 

0p  =  inpact  angle  in  degrees, 

d  =  particle  diaaeter  in  aicrons. 

P 

OONSIE 

Coefficient  in  eapirical  erosion  equation 

ail/sec 

2.22 

EXPV 

Particle  velocity  exponent  in  eapirical 

equation 

none 

6.4 

EXPRHV 

Particle  aass  flux  exponent  in  eapirical 

equation 

none 

1.0 

EIPTH 

Iapact  angle  exponent  in  eapirical 

equation 

none 

2.5 

EXPD1 

Diaaeter  first  exponent  in  empirical 

equation 

none 

1.0 

EXPD2 

Diaaeter  second  exponent  in  eapirical 

equation 

none 

6.6 

General 

Innut: 

TUB 

Tiae  period  for  calculating  total  eroaion 

depth 

sec 

- 

MODEL 

■  0  to  bypass  debris-layer  aodel 

none 

0 

calculations 
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ISAVE 


XRHOP 


Unit  on  which  particle  impact  conditions 
are  to  be  saved.  See  description  under 
IREAD  in  Section  2.2  for  parameters  saved. 
(Set  to  zero  if  conditions  are  not  to  be 
saved.) 

Particle  density  correction  factor  if  time 
particle  size-vs-mass  distribution  is  not 
used  (e.g.f  a  single  particle  size  run 
using  a  size  other  than  the  mass  mean) 


PIEROS  STAND-ALONE  CODE 


First  Card:  Title  information  containing  up  to  80  characters 


2.2.1 


DATAIN  Namelist 


Description 


Units 


1 

$DATAIN 

RT 

Nozzle  throat  radius 

inches 

- 

ZWMAX 

Ratio  of  maximum  nozzle  length/RT 

none 

- 

1 

1 

SMP 

Particle  material  density 

lbm/ft3 

- 

f 

BPS 

Particle  enthalpy  at  solidification 

n.2,  2 

ft  / sec 

- 

,  r 

ZW 

Ratio  of  nozzle  length  at  impact 

none 

- 

i 

point/RT 

f 

-f 

UG 

Gas  axial  velocity  at  impact  point 

ft/sec 

- 

, 

VG 

Gas  radial  velocity  at  impact  point 

ft/sec 

- 

; 

RHG 

Gas  density  at  impact  point 

lbm/ft3 

- 

'  *  * 

JMAX 

Total  number  of  particle  sizes  present 

none 

1 

*  4 

1 

l  “ 

at  impact  point 

DP(J) 

Particle  diameter  at  impact  point 

microns 

- 

- 

<J  «  1  to  JMAX) 

tj 

* 

UP(I) 

Particle  axial  velocity  at  impact  point 

ft/sec 

- 

C 

(J  «  1  to  JMAX) 

H- 
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VF(I)  Particle  radial  velocity  at  impact  point 
(J  -  1  to  JMAX) 

HP( J)  Particle  enthalpy  at  impact  point 

(J  -  1  to  JMAX) 

RHP(J)  Particle  density  at  impact  point 
(J  -  1  to  JMAX) 

1READ  Unit  number  from  which  previously  saved 

data  is  to  be  read.  Data  includes  aJl 
of  the  above  parameters  from  RT  to  RHP(J). 

(Set  to  zero  if  all  data  is  contained  in 
namelist.) 

$END 

2.2.2  SERODIM  Namelist 

$ER0DIN 

Same  parameters  as  $ERODIN  namelist  described  in  Section  2.1.2. 

$END 


ft/sec 
ft2/ sec2 
lbm/ft3 

none  0 


Data  for  additional  impact  points  may  be  input  in  successive  SDATAIN 
namelists,  but  the  SERODIN  namelist  is  not  repeated.  For  example,  date  for 
three  impact  points  would  be  input  in  the  following  sequence: 

TITLE 

SDATAIN 

Data  for  first  impact  point 

Send 

Serodin 

Erosion  model  parameters 

Send 

Sdatain 

Send 


Data  for  second  impact  point 


-  =« 


■*'  -  «*fc  5^>*V  AW-  'fe*  >  "‘*'  ■>“  '*K'!(' vV3 


~g*>SV3?*^ » *- -^-n<-  *~i>*-J<»  *-_.,*>«&  *  ».»**&  ??  !j<i»;^f«  '&'>  *^*»t*W 


$DATAIN 

Data  for  third  impact  point 

$END 

3.0  OUTPUT  DESCRIPTION 


3.1  ISPP  OPTION 


3.1.1  Empirical  Erosion  Model  Results 

With  the  ISPP  option,  the  nozzle  two-phase  flow  field  is  calculated, 
and  whenever  particle  impact  is  detected  at  the  nozzle  wall,  the  resulting 
erosion  is  immediately  calculated  with  the  empirical  erosion  model.  The 
erosion  results  are  then  printed  before  the  nozzle  solution  continues,  and  a 
typical  output  is  shown  in  Figure  C-4. 

At  the  initial  impact  point  for  a  particular  particle  size  group,  the 
first  line  printed  begins  "PARTICLE  IMPACT,...",  followed  by  the  wall 
coordinates  (nondimensionalized  by  the  throat  radius)  and  nozzle  area  ratio  at. 
that  point.  This  message  is  then  followed  by  a  list  of  the  mass-averaged 
properties  of  the  impinging  particles;  i.e.,  the  particle  velocity,  mass  flux, 
diameter,  and  impact  angle.  A  mass-averaged  property,  s',  for  N  number  of 
different  particle  size  groups  is  defined  as: 

N 

*$(p  V  )  a 
_  1  p  n  n 

a  *  n«l _ 

N 

:><ppyn 

n“l 


where  <ppVp)n 


is  the  particle  mass  flux  for  particle  size  group  n. 


Following  the  particle  mass-averaged  properties,  the  results  of  the 
empirical  erosion  model  in  terms  of  the  calculated  erosion  rate,  the  maximum 
erosion  depth,  and  the  axial  location  of  the  minimum  depth  are  printed.  As 
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Figure  C-4:  ISPP  Output  at  Initial  Impact  Point 
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shovra  in  Figure  C-5,  a  triangular  erosion  pattern  is  assumed,  which  is  defined 
by  the  maximum  erosion  depth  and  the  impact  angle  relative  to  the  nozzle  wall 
surface.  The  maximum  erosion  depth  is  simply  the  calculated  erosion  rate 
multiplied  by  the  exposure  time  input  in  the  namelist  $ERODIN.  If  the  erosion 
pattern  reaches  the  nozzle  exit  plane  before  the  maximum  depth  is  reached,  the 
depth  at  the  exit  plane  Yg  is  printed  instead  of  the  maximum  depth  Y^.  The 
erosion  depth  at  the  exit  plane  is  defined  as: 


where  Z^,  Z^,  and  are  the  axial  locations  of  the  impact  point,  the  nozzle 
exit  plane,  and  the  maximum  erosion  depth,  respectively. 

Downstream  of  the  initial  impact  point  for  any  one  particle  size 
group,  the  erosion  calculation  is  repeated  at  each  wall  point  in  the  nozzle 
solution.  The  typical  printout  in  Figure  C-6  shows  the  particle  mass-averaged 
properties  and  calculated  erosion  follow  the  thrust  and  mass  flow  decrements 
that  are  integrated  from  the  initial  impact  of  the  first  particle  size  group  to 
hit  the  wall. 


When  a  different  particle  size  group  first  impacts  the  wall,  the 

"PARTICLE  IMPACT . "  printout  is  repeated  prior  to  the  listing  of  the 

tnass-a\eraged  properties  and  calculated  erosion. 


At  each  wall  point  where  the  mass-averaged  properties  of  the 
impinging  particles  are  printed,  the  properties  are  stored  internally  as  a 
function  of  the  wall  point  axial  coordinate  for  subsequent  use  in  the  debris 
layer  erosion  model.  Since  this  erosion  model  requires  the  axial  distribution 
if  the  impacting  particle  properties,  instead  of  point  values  as  used  by  the 
empirical  model,  calculations  with  the  debris  layer  model  begin  only  after  the 
nozzle  solution  is  completed. 

At  the  conclusion  of  the  nozzle  flow  solution,  a  summary  table  of 
erosion  depth  versus  the  nozzle  axial  cordinate  Z/Rt>  as  shown  in  Figure  C-7, 
is  printed  following  the  usual  ISPP  "Particle  Impingement  Suanary".  At  each 
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Figure  C-6:  ISPP  Output  Downstream  of  Initial 
Impact  Point 
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Figure  C-7 :  Summary  Output  for  Erosion  Profile  Calculated 
with  Empirical  Model 
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axial  location  in  the  "Particle  Impact  Erosion  Summary",  a  check  is  made  to 
determine  the  maximum  erosion  depth  resulting  from  any  overlapping  triangular 
erosion  patterns.  Such  overlapping  can  occur  if  the  impact  angle  and/or 
maximum  erosion  depth  associated  with  each  triangular  pattern  decreases  with 


axial  distance. 


3.1.2 


Debris  Laver  Model  Results 


If  KMODEL  in  the  $ERODIN  namelist  is  greater  than  zero,  the 
calculation  of  nozzle  wall  ferosion  with  the  debris  layer  model  begins  right 
after  the  summary  of  the  empirical  erosion  model  results  is  printed.  As  shown 
in  Figure  C-8,  the  first  page  of  the  debris  layer  model  output  contains  the 
input  parameters  from  the  $ER0DIN  namelist,  followed  by  the  axial  distribution 
of  mass-averaged  particle  properties  at  the  nozzle  wall  impact  points.  This 
latter  table  is  also  saved  in  unit  ISAVE,  if  ISAVE  was  given  a  value  greater 
than  zero  in  $ERODIH .  The  axial  lengths  in  this  table,  however,  are  in  inches 
and  are  not  measured  from  the  nozzle  throat,  but  from  the  start  of  the  debris 
layer  calculation.  This  initial  position  is  either  the  location  of  the  first 
particle  impact  point,  or  if  ERMIN  is  greater  than  zero  in  $ERODIN,  the 
location  where  the  empirical  erosion  model  calculated  the  erosion  rate  to  be 
equal  to  ERMIN. 


Following  the  input  information,  the  debris  layer  model  results  are 
presented  as  a  function  of  axial  length.  The  second,  third,  fourth,  and  sixth 
columns  in  Figure  C-8  contain  calculated  debris  layer  parameters.  The  second 
column  lists  the  percentage  of  particles  approaching  the  nozzle  wall  that 
penetrate  the  debris  layer  and  impact  the  wall.  The  sixth  column  shows  the 
average  number  of  collisions  experienced  by  each  approaching  particle  within 
the  debris  layer.  Columns  three  and  four  contain  the  average  axial  velocity 
and  density  of  the  particles  that  form  the  debris  layer.  The  calculated 
erosion  rate  due  to  the  particles  penetrating  the  debris  layer  is  then 
presented  in  the  fifth  column. 


Using  the  axial  distribution  of  calculated  erosion  rates,  and  the 
assumed  erosion  pattern  (Figure  C-5),  a  summary  table  of  erosion  depth  versus 
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Figure  C-8:  Debris  Layer  Erosion  Model  Output 
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axial  distance  is  then  prepared  and  printed.  In  this  susmary,  the  axial 
distance  is  the  same  (Z/Rfc)  as  that  used  in  the  empirical  model  summary  table 
to  allow  direct  comparison  of  results. 

3.2  PIEROS  STAND-ALONE  CODE 

3.2.1  Empirical  Erosion  Hodel  Results 

With  the  PIEROS  stand-alone  erosion  prediction  code,  the 
mass-averaged  particle  properties,  the  resulting  erosion  rate  and  depth,  and 
the  axial  location  of  both  the  impact  point  and  the  maximum  erosion  depth  are 
printed  for  each  impact  location.  The  impacting  particle  properties  are  input 
either  in  the  $DATAIN  namelist  or  read  from  unit  IREAD.  The  output  format,  as 
shown  in  Figure  C-9,  is  the  same  as  that  used  within  the  ISPP  nozzle  flow 
solution.  A  summary  table  of  erosion  depth  versus  Z/R^  follows,  which  is  also 
identical  in  form  to  that  in  the  ISPP. 

3.2.2  Debris  Laver  Hodel  Results 

The  output  from  the  debris  layer  erosion  model  contained  within  the 
PIEROS  code  is  identical  to  the  debris  layer  results  obtained  with  ISPP  (Figure 
C-8) . 


3.3  ERROR  MESSAGES  AND  DIAGNOSTICS 


For  description  of  the  error  message  and  diagnostics  that  can  be 
obtained  from  the  ISPP,  the  user  is  referred  to  Section  5  of  Reference  C-l. 
The  erosion  calculations  contain  only  a  few  messages  that  slight  appear  in  the 
output,  and  these  are  described  below: 


MESSAGE 


EXPLANATION  OR  COBBKCTIVE  ACTION 


ERROR  IN  NAMELIST  ERODIN 


An  input  parameter  contained  within 
the  $ERODIN  namelist  does  not  match 
a  variable  name  in  Section  2.1.2. 
Check  input  data. 
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WARNING:  LIQUID  PHASE 
PARTICLES  PRESENT 


ZERO  EROSION  RATE  CALCULATED  FOR 
XXXX  CONSECUTIVE  STEPS 


EXCEEDED  XXXX  ELEMENTS 


TABLE  ERROR  IN  LOON 


ERROR  IN  NAMELIST  DATAIN 


The  tempeature  of  at  least  one  inpinging 
particle  size  group  is  above  the  nelting 
temperature  of  the  particle  material. 

See  Section  4  for  recommendations  for 
these  conditions. 

The  debris  layer  erosion  model 
calculated  erosion  rates  (less  than  0.01 
mil/sec)  for  XXXX  consecutive  debris 
layer  elements  and  the  run  was 
terminated.  This  prevents  excessive 
computer  time  usage  and  is  controlled  by 
the  input  parameter  IEMAX  in  the  $ER0DIN 
name  list . 

The  maximum  number  of  debris  layer 
elements  used  in  the  debris  layer 
erosion  calculation  was  exceeded  and  the 
run  was  terminated.  This  prevents 
excessive  computer  time  usage  and  is 
controlled  by  the  input  parameter  MAXY 
in  the  $ER0DIN  namelist. 

An  axial  location  vas  not  contained 
vithin  the  input  table  of  particle 
conditions  versus  length  and  the  run 
vas  terminated.  Check  input  data. 

An  input  parameter  contained  vithin  the 
$DATAIN  namelist  for  the  PIEROS  code 
does  not  match  a  variable  name  in 
Section  2.2.1.  Check  input  data. 
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4.0 


RECOMMENDED  PRACTICES  AND  GOIDELIKES 


4.1 


ISPP 


For  general  recommended  practices  and  guidelines  £or  using  ISPP,  the 
user  is  referred  to  Section  4  of  Reference  C*-l.  A  few  additional 
recommendations,  however,  can  be  made  to  help  give  more  reasonable  particle 
erosion  resuts. 


4.1.1 


Particle  Size-vs-Mass  Distribution 


It  is  strongly  recommended  that  the  low  end  of  the  particle 
size-vs-mass  distribution  be  well  represented  since  the  majority  of  the 
particle  impact  will  be  produced  by  the  small  sizes.  The  ISPP  has  been 
modified  to  allow  this  distribution  to  be  specified  by  inputting  mass  fractions 
for  which  corresponding  particle  sizes  will  be  calculated  internally.  For 
cases  run  to  date,  the  sir.  particle  size  distribution  corresponding  to  WPWT  * 
.01,  .02,  .04,  .08,  .15,  and  .70  has  proven  to  be  successful.  In  larger 
nozzles,  some  problems  in  the  nozzle  flow  solution  may  be  experienced  with  the 
smallest  particles.  In  such  cases,  the  distribution  can  be  modified  slightly 
by  using  five  sizes  corresponding  to  WPWT  *  .02,  .04,  .08,  .16,  and  .70,  or 
four  sizes  corresponding  to  WPWT  *  .04,  .08,  .16,  and  .72.  Mote  that  the  sum 
of  the  WPWT  values  must  equal  1.0. 


The  ISPP  computer  run  time  increases  almost  linearly  with  the  number 
of  particle  size  groups  used.  Therefore,  in  some  cases  it  may  be  worthwhile  to 
make  a  single  particle  ISPP  run  with  a  small  particle  size  representative  of 
the  lower  end  of  the  size-vs-mass  distribution  to  determine  if  particle  impact 
will  occur  in  a  given  nozzle.  Single  particle  results  have  correlated  fairly 
well  with  those  for  multi-particle  runs  as  far  as  impact  location  is  concerned. 
However,  to  obtain  reasonable  calculated  erosion  results,  it  is  necessary  that 
the  proper  mass  fraction  be  assigned  to  the  single  particle  size  by  using  the 
variable  XRH0P  (see  Section  2.1.2). 
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4.1.2 


-*  rt 


TD2P  Control  Parameters 


Two  control  parameters  input  within  the  $TD2  namelist  that  are 
important  to  the  particle  trajectory,  and  thus  impingement,  calculations,  are 
DHII  and  DCR.  DH1T  specifies  the  maximum  separation  distance  that  is  allowed 
between  a  particle  streamline  and  nozzle  wall  before  impact  is  assumed  to 
occur.  Too  large  a  value  for  DHIT  forces  an  artificial  particle  impact  at  an 
axial  location  too  far  upstream.  On  the  other  hand,  too  small  a  value  can  make 
the  impact  location  calculation  too  dependent  on  small  numerical  errors  in  the 
flow  field  solution,  particularly  at  low  impact  angles.  DCR  specifies  the 
maximum  separation  distance  that  is  allowed  between  adjacent  particle 
streamlines  for  different  particle  size  groups  before  the  streamlines  are 
assumed  to  cross.  Too  large  a  value  for  DCR  can  artificially  shift  layer 
particles  closer  to  the  wall,  greatly  changing  any  subsequent  impact 
conditions.  Like  DHIT,  too  small  a  value  for  DCR  can  lead  to  numerical 
difficulties.  Based  on  past  experience  with  numerous  TD2P  calculations,  it  is 
recouuended  that  the  values  of  0.005  and  0.010  be  used  for  DHIT  and  DCR, 
respectively,  for  nozzles  with  throat  radii  in  the  1.5  to  6.0  inch  range.  For 
larger  nozzles,  these  parameters,  which  are  actually  separation  distances 
divided  by  the  throat  radius,  should  be  reduced  to  maintain  approximately  the 
came  physical  separation  distances. 

Other  TD2P  control  parameters  that  can  affect  the  calculated  particle 
trajectories  are  NILP,  the  number  of  points  along  the  HOC  start line,  and  both 
DR  and  DL,  the  maximum  width  between  solution  points  along  a  right-running  and 
left-running  characteristic,  respectively.  Too  small  a  value  for  NILP  and  too 
large  a  value  for  DR  and  DL  can  lead  to  solution  inaccuracies  and  possible 
failure,  while  vary  large  computer  run  times  can  result  with  a  large  NILP  and 
small  DR  and  DL.  It  was  found  in  studies  conducted  during  this  program  that 
the  default  values  of  NILP  «  40  and  DR  *  DL  ■  0.1  gave  the  best  results,  and 
are  therefore  recomended. 


-a 


I 


Be 


4.1.3  Transonic  Solution 


would  definitely  be  preferred  since  the  ATA  is  limited  to  conical  inlet  shapes. 
The  FCT,  however,  is  a  more  complex  solution  requiring  greater  amounts  of  input 
preparation  and  computer  run  time,  and  leading  to  more  potential  problems  in 
achieving  a  successful  solution.  If  the  ATA  i6  used,  the  nozzle  conical  inlet 
angle  and  upstream  wall  radius  of  curvature  should  be  selected  to  best 
represent  the  actual  nozzle  inlet  geometry.  For  nonsubmerged  nozzles,  inlet 
angles  of  30°  or  less  should  be  used  to  give  representative  particle 
trajectories,  and  for  submerged  nozzles,  a  60°  inlet  angle  is  recommended. 

4.2  EROSION  CODES 


4.2.1 


Use  of  Empirical  and  Debris  Laver  Erosion  Models 


The  debris  layer  erosion  model  can  use  significant  amounts  of 
computer  time  for  cases  where  relatively  large  number  of  collisions  occur 
within  the  debris  layer.  It  is  therefore  recommended  when  using  the  ISPP 
option  that  the  debris  layer  calculations  not  be  made  in  the  seme  run  as  the 
nozzle  flow  solution  by  setting  KMODEL*0  in  the  $ERODIN  namelist,  and  that  the 
particle  impact  conditions  be  saved  by  setting  ISAVE  in  $ER0DIN  to  the  unit 
number  on  which  the  data  is  to  be  written.  In  this  manner,  the  results  of  the 
empirical  erosion  model  can  be  examined  to  determine  if  debris  layer 
calculations  are  necessary.  If  so  desired,  the  debris  layer  model  can  then  be 
run  with  the  PIEROS  code  and  the  data  from  unit  ISAVE. 

In  evaluating  the  results  from  the  empirical  erosion  model,  the 
calculated  erosion  rates  would  be  expected  to  be  valid  if  the  particle  flow  and 
impact  conditions  fall  within  the  followiag  ranges: 


Particle  velocities  from  4000  to  11000  ft/s, 

2 

Particle  mass  fluxes  less  than  20  lbm/s*ft  , 

Particle  diameters  less  than  10  pm. 

Particle  temperatures  less  than  the  particle  material  melting 
temperature,  and 


Impact  angles  less  than  10  . 
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For  particle  flow  and  impact  conditions  outside  the  above  ranges,  the 

debris  layer  model  should  be  used,  unless  the  particle  mass  flux  is  below  1.0 
2 

lbm/s’ft  .  If  the  particles  are  in  the  liquid  phase,  a  first  order  correction 
to  the  erosion  rates  calculated  by  the  debris  layer  model  can  be  made  by 
dividing  the  rates  by  a  factor  of  20. 

4.2.2  Nozzle  Wall  Material 

The  erosion  models  were  developed  primarily  for  carbon-carbon  and 
graphite  nozzle  wall  materials.  For  the  empirical  model,  a  first-order 
correction  for  carbon-phenolic  materials  based  on  material  density  has  been 
derived,  but  this  correction  would  not  be  expected  to  be  valid  for  other 
materials.  The  debris  layer  model  does  use  the  yield  shear  stress  and  density 
of  the  target  material  (ST  and  RHOTAR  in  the  $ER0D1N  namelist),  and  this 
theoretically  could  give  representative  erosion  results  for  other  materials. 
However,  the  debris  layer  model  results  have  only  been  evaluated  for  the 
carbon-carbon  and  graphite  materials  to  date. 

4.2.3  Scale  Effects 

Correlation  studies  have  shown  that  the  erosion  models  gave  the  best 
rsults  for  nozzles  with  throat  radii  in  the  range  from  1.5  to  6.0  inches.  For 
smaller  nozzles,  the  erosion  was  usually  overpredicted,  and  for  larger  nozzles 
underpredicted.  This  effect  is  thought  to  be  related  more  to  the  two-phase 
flow  solution  rather  than  the  erosion  models  themselves. 

4.2.4  Erosion  Model  Input  Parameters 

As  shown  in  Section  2.1.2,  default  values  based  on  previous 
correlations  are  contained  within  the  erosion  codes  for  most  of  the  input 
parameters.  The  only  parameter  that  always  needs  to  be  input  is  TIME,  the  time 
period  used  to  convert  the  calculated  erosion  rates  to  erosion  depths. 
However,  the  default  values  can  be  overridden  if  different  options  are  designed 
(such  as  saving  data  on  unit  ISAVE),  or  if  the  correlation  with  a  specific  set 
of  erosion  data  can  be  improved  by  adjusting  one  or  more  of  the  parameters  for 
either  of  the  erosion  models. 


5.0 


SAMPLE  CASES 


Three  sample  cases  are  presented  to  illustrate  the  use  of  the 
particle  impact  erosion  models.  The  first  is  a  complete  ISPP  r<_n  with  the 
complete  nozzle  flow  solution  and  the  erosion  calculated  at  each  wall  point 
where  impingement  occurs.  The  particle  flow  conditions  at  the  impact  points 
were  saved  for  later  input  to  the  stand-alone  code  PIEROS,  which  is  illustrated 
in  the  second  sample  case.  The  third  sample  case  also  used  PIEROS,  but  with 
all  the  particle  conditions  input  through  the  namelist  $DATAIN. 

For  all  three  cases,  the  SD001  motor  and  nozzle  configuration  was 
used  since  it  experienced  severe  particle  erosion  in  the  nozzle  exit  cone. 

5.1  CASE  1:  ISPP  OPTION  WITH  EROSION  MODELS 

Sample  Case  <  includes  a  full  two-phase  nozzle  flow  solution  with 
particle  impact  erosion  calculated  with  both  the  empirical  and  debris  lcyer 
model  within  the  same  run. 

The  control  cards  used  for  this  run  on  a  UNIVAC  1108  computer  are 
shown  in  Table  C-l.  In  this  case,  the  gas  and  particle  properties  were  input 
directly;  but  if  the  properties  are  to  be  calculated  with  the  ODE  code, 
additional  control  cards  are  required  to  assign  a  thermodynamic  property  data 
file  and  equate  it  to  unit  25  (see  Section  2.10  of  Reference  C-l).  At  the 
completion  of  the  ISPP  run,  the  impacting  particle  conditions  were  written  on 
unit  10  and  saved  as  an  element  in  a  user  file  with  the  @ELT  and  @ADD  cards. 
Note  that  to  create  this  element,  the  ISPP  run  cannot  terminate  on  an  error 
message.  This  was  the  reason  the  main  program  DRIVER  has  to  be  slightly 
modified. 


The  input  data  for  Case  1  is  listed  in  Table  C-2.  Note  that  in  the 
$SPP  namelist,  only  the  TD2P  flag  is  nonzero,  and  only  the  nozzle  geometry  data 
is  included.  In  the  $TD2  namelist,  all  the  required  gas  and  particle 
properties  are  input  in  the  first  four  lines.  If  the  ODE  were  to  be  used  to 
calculate  these  properties,  the  ODE  flag  in  the  $SPP  namelist  would  be  set  to 
1,  and  additional  input  included  before  the  $TD2  namelist.  A  third  option  can 
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@XQT 
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■ 
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S 

i 

§ 

SPP. 

10. 

SPP.SPPA 


Assign  program  file  for  ISPP 

Assign  temperory  file  for  unit  10  (i.e.  ISAVE) 

Execute  absolute  element  SPPA  in  program  file  SPP. 


DATA 


@ELT, IL  SPP.DATA10 
@ADD,P  10. 

@FIN 

@@ 


Create  element  DATA10  in  program 

file  SPP.  (or  any  other  file)  with  data  written 

on  unit  10  (i.e.  ISAVE)  during  program  execution 


Table  C-l:  UNIVAC  1108  Control  Cards  for  Sample  Case  1 
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3^ 


ssSaSs. 


TITLE1  SD001 
SPP 

sspp; 

OOE=0.*TD2P=1.,TBL=0., 

INLET=0,RWTU=2.16,THETAI=30*, 

IWALL=4,RWTD=1.33,THETA=S0.B,THE=12.5» 

RS<2)=7, 315, 11. 214, 15. 942, 20. 167, 24, 726, 28. 137, 31. 945, 
ZS(2)=7. 490, 15. 602, 27. 366, 40. 296, 56. 713, 70. 452, 87. 02* 
NWS=8,«ZN0RM=3.B81, 

RSI=3.B81, 

SEND 

*TD2 

TG0=6345.,GMG0=. 6342-4, CAPN=. 67756, CPL=11288., 
TPM=4167.,EPM=6. 3752, UGH=7951,,DHH=1. 25+7, PR=. 4617, 
CPS0L=»149.7, 

WPWST=*45762, 

0L=.1,0R=.1,THIW=1., 

NIL®=40, 

NPG=6, WpwT=»Ol, •02,«04,o06,»l5,*70, 

PC-BOO., 

DHIT=.005, 

DCR=. 01, 

SEND 
serdoin 
TIM£i=4*  , 

ISAtfE=10, 

SEND 


Table  C-2:  Input  Data  for  Sample  Case  1 
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be  called  by  setting  ODE  equal  to  2  end  reading  in  the  property  data  from  a 
file  created  during  a  previous  ODE  run.  For  details  on  using  the  ODE  *  1  and 
ODE  »  2  options,  the  user  is  referred  to  Reference  C-l. 

The  $TD2  namelist  in  Table  C-2  also  includes  the  recommended  values 
for  the  control  parameters  DL,  DR,  NIL?,  DHIT,  and  DCR,  and  the  recommended 
particle  mass  fractions  to  be  used  for  calculating  the  particle  size 
distribution.  The  ATA  transonic  solution  for  conical  inlets  was  used  for  this 
case.  However,  if  the  FCT  transonic  solution  was  desired,  additional  input 
would  be  included  in  the  $TD2  namelist  (see  Reference  C-l). 

In  this  case,  most  of  the  default  values  for  the  erosion  model  input 
parameters  were  used,  and  thus  only  the  erosion  time  period  and  the  ISAVE  unit 
number  for  serving  the  calculated  particle  conditions  at  impact  are  input  in 
the  namelist  $ER0D1N. 

Selected  portions  of  the  output  for  Case  1  are  presented  in  Table 
C-3.  The  first  21  pages  of  this  table  show  the  standard  1SPP  output,  including 
the  input  from  the  SPP  namelist,  the  gas  and  particle  properties,  the 
calculated  nozzle  contour,  the  calculated  particle  size-vs-mass  distribution, 
the  HOC  start  line  determined  from  the  ATA  transonic  solution,  and  some  of  the 
HOC  supersonic  two-phase  flow  solution.  Page  57  show;’  the  output  where  the 
first  particle  impact  is  detected  for  Particle  Group  1  with  a  1.278  micron 
diameter,  and  the  resulting  erosion  calculated  with  the  empirical  erosion 
model.  The  following  page  shows  the  erosion  output  for  the  same  particle  size 
at  a  subsequent  wall  point.  The  results  for  the  first  impact  point  for 
Particle  Group  2  (1.782  micron  diameter)  are  then  shown  on  the  following  page 
(page  59).  The  calculated  erosion  at  this  point  is  based  on  the  mass-averaged 
properties  of  the  combined  particle  Groups  1  and  2.  On  subsequent  pages,  the 
appearance  of  the  words  'PARTICLE  IMPACT,  HALL  (R,Z)  ■  ....,'  indicates  the 
initial  impact  point  for  a  different  particle  size  group. 

The  nozzle  flow  solution  for  Case  1  is  concluded  on  page  77,  followed 
by  a  PARTICLE  IMPINGEMENT  SUMMARY  that  shows  the  net  mass  flow  and  thrust 
decrements  associated  with  each  particle  size  group  that  imputed  the  wall 


within  the  nozzle.  The  calculated  erosion  depths  iron  the  empirical  model  are 
then  summarized  as  a  function  of  the  nondimensional  axial  distance  from  the 
throat,  Z/Rfc. 

The  input  data  for  the  debris  layer  erosion  model  then  follows  on 
page  79.  The  table  of  impacting  particle  conditions  versus  length  was  stored 
internally  and  also  written  on  unit  10,  corresponding  to  the  ISAVE  value  in  the 
$ER0D1N  namelist.  The  debris  layer  results  are  shown  on  page  80,  the  last  page 
of  Table  C-3.  In  this  case,  the  density  of  the  particles  in  the  calculated 
debris  layer  was  large  enough  to  protect  the  nozzle  wall  downstream  of  the 
second  print  location,  resulting  in  the  calculation  being  terminated  when  a 
zero  (actually  less  than  0.01  mil/src)  erosion  rate  was  calculated  for  100 
consecutive  steps.  Termination  for  more  or  less  consecutive  steps  can  be 
controlled  by  inputting  the  desired  maximum  number  of  steps  as  IEMAX  in  the 
$ERODIN  namelist. 

The  last  portion  of  the  Case  1  output  is  a  suimnary  of  erosion  depths 
calculated  with  the  debris  layer  model  as  a  function  of  Z/R  .  Comparing  this 

l 

table  with  the  corresponding  table  containing  the  empirical  erosion  model 
results  on  page  78,  it  can  be  seen  that  the  debris  layer  model  apparently 
overpredicts  the  shielding  effect  for  this  case,  since  the  calculated  erosion 
is  much  less  than  the  empirical  model  results,  which  correlate  well  with  data. 
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5.2  CASE  2:  PIEROS  WITH  DATA  FROM  ISPP 

The  second  sample  used  illustrates  the  use  of  the  PIEROS  stand-alone 
erosion  code  using  the  particle  conditions  saved  on  unit  10  during  the 
execution  of  Case  1.  This  option  is  useful  if  the  run  time  of  the  ISPP  run  is 
to  be  minimized  by  not  including  the  debris  lfver  erosion  model  calculations  in 
the  same  run.  Also,  input  parameters  for  the  empirical  eros  n  model  and/or 
the  debris  layer  model  can  be  changed  without  re-running  the  ISPP  nozzle  flew 
solution.  As  an  example,  the  erosion  characteristics  of  different  candidate 
nozzle  wall  materials  could  be  evaluated  by  changing  the  target  material 
properties  in  subsequent  PIEROS  runs  using  the  same  input  impacting  particle 
conditions . 


The  UNIVAC  1108  control  cards  for  Case  2  are  listed  'n  Table  C-4. 
Here  the  conversion  of  the  saved  particle  conditions  from  an  element  in  a  users 
file  to  the  unit  to  V 2  read  (specified  by  IREAD  in  the  $DATAIN  namelist)  is 
accomplished  with  the  @DATA  and  6ADD  control  cards. 

The  input  for  Case  2  is  shown  in  Table  C-5  to  be  very  simple, 
consisting  only  of  the  IREAD  description  in  the  $DATAIN  namelist,  and  the 
erosion  time  period  in  $ERODIN.  Again,  as  in  Case  1,  the  default  values  for 
most  of  the  erosion  model  input  are  being  utilized.  However,  in  subsequent 
runs,  these  default  values  could  be  overwritten  by  including  new  values  in  the 
namelists. 


The  Case  2  output  is  presented  in  Table  C-6,  covering  pages  84  to  93. 
The  first  three  pages  show  the  particle  conditions  saved  from  Case  1  (and 
described  in  Section  3.  ).  On  page  87,  the  erosion  calculated  with  the 

empirical  erosion  model  is  given  for  each  input  point  using  the  same  output 
format  as  Case  1,  where  the  erosion  results  are  printed  within  the  nozzle  flow 
solution.  The  summary  of  the  erosion  depths  from  the  empirical  model  and  the 
debris  layer  erosion  model  results  follow  in  the  same  order  and  format  as  Case 


"W^  -ta  '-'V^.''-*  ^‘Sh^’wSS* 


@ASG,AX  PIE. 

@ASG,AX  SPP. 

@ASG,T  10. 

@DATA,IL  10. 

@ ADD, PD  SPP.DATA10 

@XQT  PIE. PIEROS 

DATA 


Assign  program  file  for  PIEROS  stand-alone  :ode 
Assign  program  file  containing  data  saved  from 
previous  I SPP  run. 

Assign  temporary  file  for  unit  10  (i.e.  IREAD) 
Assign  previously  saved  data  file  as 
Unit  10  (i.e.  IREAD) 

Execute  absolute  element  PIEROS  in  program 
file  PIE. 


@F1N 

@@  Remote  job  entry  end  of  file 


Table  C-4:  UNIVAC  1108  Control  Cards  for  Sample  Case  2 


SDOOl  TEST  case 

SDATAIN 

IREAD=10* 

SEND; 

SERDDIN 

TIM£t=4.« 

SEND 
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MASS-AVCRAGEO  ^"OPEATICS  0*  1MPINGIN8  PARTICLES  AT  Z/RT  s  JZ.S9I  t 

VELOCITY  «  10*73.05  PT/SEC*  IASS  PUN  •  1.0*  LP VSEC*PT«*2.  PARTICLE  01AM.  x  3.29%  fllCRONS,  IPP«T  ANGLE  «  t.oc#  OES 

CALCULATED  EMPIRICAL  EROSION  RATE! 

14.0*3  HILS/SEC.  RESULTING  IN  A  MAXIMUM  EROSION  DEPTH  OP  *.19*  NILS  AT  Z/R»  «  **.■(** 
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5.3 


CASE  3:  PIEROS  STAND-ALONE  CODE 


Case  3  shows  that  the  esipirical  and  debris  layer  models  can  be  used 
conpletely  independently  of  the  ISPP.  If  the  impacting  particle  flow 
conditions  are  known  from  another  source,  or  if  sensitivity  studies  are  being 
conducted  to  determine  the  effect  of  a  particular  particle  flow  condition  on 
erosion,  the  conditions  may  be  input  directly  to  the  PIEROS  stand-alone  code 
within  the  $DATAIN  namelist. 

The  UNIVAC  1108  control  cards  and  the  input  for  Case  3  are  shown  in 
Tables  C-7  and  C-8,  respectively.  Notice  in  Table  C-8  that  if  an  axial 
distribution  of  impacting  particle  conditions  is  to  be  input,  an  additional 
§DATAIN  namelist  is  required  for  each  impact  point.  However,  the  $ERODIN 
namelist  is  required  only  once,  following  the  first  $DATAIN* 

The  output  format  for  Case  3  is  shown  in  Table  C-9  to  have  the  same 
sequence  and  format  as  the  Case  2  output  in  Table  C-6. 


@ASG,AX  PIE.  Assign  program  file  for  PIEROS 

Stand-alone  erosion  code 

;?XQT  PIE. PIEROS  Execute  absolute  element  PIEROS 

in  program  file  PIE. 

DATA 

@FIN 

@@  Remote  job  entry  end  of  file. 


Table  C-7:  UNIVAC  1108  Control  Cards  for  Sample  Case  3 

fi 
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SDODl  TEST  CASE 
sdatain. 

RT=3.B8l,ZUHAXs22.922,snP=247.7,HPS=l.* 
ZW=9,853*UG=9787.*VGs:2986.,RHS=. 00318* 
JMAX=1*UP=9636.*VP=3Q85**RHP=1.&11E-5*HP=0.*DP=1.27?* 

SENG 

serodin 

C3=5.  ,RHOTAR=120.  •  Y=200.  *ST=2200.  ,XPRINT=2.  •  TlM£=8.  ,KNDD.-  .i  t 
S1AXT=8000* 

ERMlM=l.* 

iSAVEslO* 

SEND1 

sdatain 

ZW-l0.689*UG=9859.  ,VG=2923.  *RHG=.  00131* 

JMAX=2*UP=9756,*9702.*VP=2977.*2998.*RHP=1.833E*5,2.776E-5* 
HP=2*0 • * DP=1 • 277* 1 . 782 * 

SEND 

sdatain 

ZWsll.‘»a&,UG=9968,.VG=2863.,RHG=. 00259, 

JMAX=3«UP=9856,  *9808.  «97bO.*VP=29l9.  *293*1.,  2958.* 
RHP=1.282E-5,2.535E-5*5.035E-5*0P=1.277*1.782,2.319* 

HP=0. *0. ,0. « 

SEND 

sdatain 

'  ZW=12.653*UG=10092. • V6S2785. , RHG=. 00227* 

JMAXs:8«UP=998«»#. 9936.. 9881. *9799.* VP=2836.  ,2857.  *2878. *2905* • 
RHP=1 .  U5E-5 . 2 . 200E-5  «  8 . 333E-5 . 8 . 677E-5  .  HP=8*0  .  i 
DP=1. 277, 1.782, 2. 318. 3. 028, 

SEND 

SDATAIN 

ZUsl8.969*UG=10296.«VGs2639,,RH5-.00181* 
v*HAX=5»UP=10189.  *10139.  *10085.  *10008.  *9662.* 

VP=2703. ,2725. ,2785*  *2771. *2809. *HP=5*0,  * 
RHP=a.70*E-6»1.709E-5,3.351E-5«6.606E-5*1.278E'*8, 

DPrl. 277,1. 782,2.319, 3,028*8, 088* 

SEND 


Table  C-8:  Input  Data  for  Sample  Case  3 
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